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TRANSPORT OF INFRARED RADIATION IN 
CU601DAL CLOUDS 

Harshvardhan 
James A Weinman 

Roger Davies 

abstract 

The transport of infrared radiation in a single cuboidal cloud has been modeled using a verti* 
cal rwo~stream approximation. Computations have been made at lO/xm for a Deirmendjian (1969) 
C-1 water cloud of single scattering albedo, Z » 0.638 and asymmetry parameter, g « 0.865. Results 
indicate that the emittance of the top face of the model cloud is always less than that for a plane par* 
allel cloud of the same optical depth. The hemispheric flux escaping from the cloud top has a gra- 
dient from the center to the edges which brighten when the cloud is over warmer ground. Cool- 
ing rate calculations in the 8-’13.6Mm region show that there is cooling out of the sides of the 
cloud at all levels even when there is heating of the core from the ground below. 

The radiances exiting from model cuboidal clouds were computed by path integration over 
the source function obtained with the two-stream approximation. Results suggest that the 
brightness temperature measured from finite clouds will overestimate the cloud top temperature. 

Some key results of the model have been compared with .Monte Carlo simulations. Overall 
errors in flux and radiance average a few degrees for most cases. 
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TRANSPORT OF INFRARED RADIATION IN 
CUBOIDAL CLOUDS 


1. INTRODUCTION 

The role of clouds in determining the radiation field in general circulation and climate 
models has been widely considered. It has also been realized that the finite horizontal dimen* 
sions of individual clouds in convective cloud fields cannot be treated in the same manner as 
plane parallel clouds. The scattering of visible radiation by finite clouds has been studied by 
McKee and Cox (1974) and others, prim?nl/ using Monte Carlo methods, and by Davies and 
Weinman (1977), Rarkstrom and Arduim (1977) and Davies (1978) who developed analytical 
solutions to this problem. Microwave radiances from finite rain clouds were computed by Wein- 
man and Davies (1978) using an analytical method that was compared with Monte Carlo results. 
The question of infrared radiances from finite clouds has been addressed only recently by Liou 
and Ou (1979) and EUingson and Kolczynski (1980); however, a number of considerations still 
require attention. 

The purpose of this paper is to establish a simple analytical method that can be used to 
compute infrared radiation quantities in isolated finite clouds of cuboidal shape. This will be 
accomplished by comparing some key results with Monte Carlo computatic:^. The analytical 
technique is such that allowances for vertical temperature variations can be made and in the 
future, the isolated model cloud may be embedded in a participating atmosphere. We have also 
not considered the effect of nei^boting clouds. 


1 


2, VERTICAL TWO-STREAM APPROXIMATION 


The equation of transfer in three dimensions cannot be solved analytically for arbitrary scat- 
tering. Various approximations to the equation and boundary conditions have therefore been 
suggested and utilized to obtain results for particular problems. Two recent applications in at- 
mospheric radiative transfer are those of Davies (1978) who used the Delta-Eddington approxi- 
mation in the visible and Liou and Ou (1979) who used the Eddington approximation in the 
infrared. Weinman and Davies (1978) have also computed microwave radiances emerging from rain 
clouds using the Eddington approximation. A problem with the Eddington approximation is that it 
can yield unrealistic results such as negative reflectance and values of emittance exceeding unity for 
certain ranges of the single scattering albedo, u, and asymmetry parameter, g. Lyzenga (1973> 
has shown that for a serai-infinite plane parallel layer, negative values of the diffuse reflectance 
are obtained in the Eddington approximation when 

g > (4S5 - 1) (1) 

Unfortunately the range of parameters typical of clouds satisfies the above inequality because g ^ 

0.7 to 0.8 and 0.3 to 0.6 in the infrared; hence, the Eddington approximation is inappropriate. 

The two-stream approximation avoids the pitfall that Liou and Ou (1979) encountered, yet 
it yields similar governing equations. Recognizing that stratification within a horizontally finite 
cloud element would stiU be in the vertical direction, it is possible to introduce a vertical two- 
stream approximation to the three dimet^sional i idiative transfer problem. The equation of trans- 
fer in three dimensions for a homogeneous cloud may be written in terms of the radiance, 1. as 

SI SI SI . 

v COS ^ Sin 0 M T- ^ - I + (1 “ 25) B + J, (2) 

8x Sy Sz 

The horizontal and vertical axes are designated x, y and z, where z increases in the upward 
direction (see Fig. 1). The independent variables are in dimensionless optical thickness, r, 
with the extinction coefficient, k, absorbed into the L.H.S. of the equation. The azimuth and 
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zenith angles, 0 and $, where m " cos ^ and v * $in 9, are measured with respect to the x and z 
axes, respectively. The scattering source func*^ ^n is 

(x, y* ^ y y p ^ 

c3 is the albedo for single scattering, p (m. m'. 0') is the phase function for scattering from 

direction n\ 0' into m> 0> The Planck function B may be linearized about its value at the top of 
the cloud, such that 

B - Bo + jte’ (4) 

where 0 is a constant. However, for the purpose of this study the cloud is assumed to be 
isothermal, such that B • Bq. All quantities are monochromatic and cloud properties are con- 
sidered to be uniform throughout. For illustrative purposes, a wavelength of 1 0/btni is used in 
most of this work. 


The two-stream approximation assumes that two quantities, I'*’ and r represent the radiance 
and further that they are identified with the value of the radiance at m ® which are the 

first-order gaussian quadrature points. Thus the scattering source function (3) may be written 
for the positive (upward) hemisphere as 

~ y" y I (x, y, z, n\ 0') p 0; m'. 0'j d>i' d0' 


^ f i /r \ 

Ll'‘'U.y.z.0')pl y j.0:yj.0'l 
+ r(x,y,z,0')p ^y^.0;- y^,0'^ j d0' 


(5) 


where I* is identified with m ■ ^ • If the phase function is ex- 

panded in a series consisting of the zero and first order Legendre functions, then 


p (/I, 0 :m', 0') * 1 ♦ 3g [mm' + w' cos 10 - 0')) 


( 6 ) 


- 


where g is the asymmetry parameter. Substituting for ju and n', the R.K.S. of (5) reduces to 
I g 2gcos(0 - ^')1 + r(l - g + 2gcos(0 - 0^)]} d0* 


Similarly for the downward hemisphere, we have 

•2ff 


7” / {l“ (1 + g + 2g COS(0 - 0')] + - g + 2gcos(0 - 0')l} d0' 

The equation of transfer in the vertical two-stream approximation follows from (2) 

/T dl+ /T dl+ . 31+ H. -txo 

y? IT * /t -sr * V T 17 * ■ ‘^ * '■ * ® 

25 

+ — j {l+(l + g 2gcos(0 - 0')] 


+ r[l - g + 2gcos(0 - 0')) } d0' 


and 


/Tar /i 3r /T . ar ^ 

■ >/Tir w T 17 17 ■ ■ 


a 

+ — / { r [ l + g + 2g cos (0 - 0')1 


+ I+tl - g + 2gcos(0 - 0'))} d0' 


(7a) 


(7b) 
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3. GOVERNING EQUATIONS AND BOUNDARY CONDITIONS 

Eqs. (7a) and (7b) may be reduced by an expansion which employs the zero and first 
order Legendre functions: 

I (x, y, z, M. 0 ) ■ Ic (X, y. z) + li (X, y, z) M + Ix (x, y. zW cos ^ (x, y, z) v sin 0 

In the two-stream approximation we will have 

and 


setting 


r(x,y,2,0) • 


i; 


-yf 


i; + 


Iicos« +yjlyjin# 


'x-yf'i 


yields 


and 


Further, letting 


and 


I'*’ ■ I* + I^ cos ^ -r ly sin 


r » + I^ cos ^ + ly sin 


‘o - J (? ♦ 1? 


I, - - (I* - 1,-) 


(8a) 

(8b) 


the scattering source function in (7a) and (7b) may be written as 
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cSIo *•’ c3gli * &}gdxCOS0 -t* IySin0) positive hemisphere 
cSIq - &}gl| c3g(lxcos0 '•* lySin^) negative hemisphere 

I 

Now, integration of Eqs. (7a) and (7b) by r— I dd after substitution of (8a) and (8b) 

2ff Jq 

yields the following two equations 


(9a) 


and 


. /LE* /Lf2k.?kV-ii 

V 3 V 6 \ 3x by J * 


+ (1 - as)Bo + aio - sjgii 


(9b) 


Adding (9a) and (9b) yields 


v/y *'A{^ 


(10) 


It can be shown that (10) satisfies the flux integral relation for this problem if flrst-order gaussian 
quadrature is used to evaluate integrals over n. This equation is found in Chandrasekhar, ( 1960). 




(ID 


where Fnei is the net flux vector. Subtracting (9b) from (9a) yields 
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Integration of (7a) and (7b) by 




COS d d0 yields 


and 


-L 

> Jo 

1 A!k.i /!?£.. ii.sgii 

2V3dz 2V3dx 2 2 

.i /Lilx.i 

2 V 3 3z 2 V 3 3x 2 * 2 


( 12 ) 


(13a) 


(I3b) 
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Adding » 13 j) and (1 3b I w< have 


Vylf *M1 -ag)i, (14) 

Similarly it may be shown that 

yf^--o-3i>v (•« 

Combining (10), (U), (12). (14) and (15) yields 


7^ lo - - Bo) 06) 

where 

X ■ i3(l - c3)(l - SSg)]^^ 

is the inverse diffusion length. The boundary conditions for (16) are obtained from the flux 
conditions at the face of the cloud. Again, using the two-stream technique the fluxes may be 
written as 


Fi-.z- 
^x+, z+ 
z- 

^X+. X* 


If 

•'o •'0 
*2>r yO 


I**'MdMdd * ffdo + Ii) ■ ff 


J PMdfidd * -ffdo - Id • ho jT ~ ) 

ff 

^ J I"*" V cos 0 dM dd ■ “ ^ 


J* r cos 0 d/i dd ■ — ^ly + Ii + — 


■ Px+. z+ + 


/ ir \ / 1 dlo\ 

Fx+ 2- ■ * ho ■** “• Ix “ lo 

x+.z- y 4 y V h 3x/ 


(17a) 
(17b) 
(17c) 
(I7d) 
(1 7c) 


Similarly 


Fx-.x- * -» 

/ i aio\ 

it) 

(170 

Fy,.y, - » 1 

/ h dyj 

(17g) 

Fy..y. • -* 

/ 1 aioN 

(yzi;) 

(17h) 


Here, h' • /T (I Of) and h - 7 <' ’ ®»)- "*««» «« iUuetnted in Figure I 

following Davies (1978). It is assumed that there is no flux incident at the top of the cloud and 
the bottom receives radiation only from the groi .1 emitting at its black body temperature, T^. 
The ground flux is and the ground reflectivity is zero. The sides of the cuboidal cloud re 
ceive only the component of the flux emitted by the ground that is normal to the face. The use 
of flux boundary conditions leads to the loss of all directional information pertaining to the in- 
coming radiation. This leads to slightly erroneous results for the radiance, a matter that is dealt 
with in Section 5 when Monte Carlo comparisons are made. Following the sign convention that 
flux is positive in the positive direction, the boundary conditions may now be written for a cloud 


(18a) 


of optical thickness z* and equal sides, s. 

These are 

(u 

1 

^ «— 

dIo\ 


▼ 

h' 

bz I 

(u 

1 

^lo\ 


'7 

bz I 

(u 

1 

bJo\ 

V 

h 

dx4 

(u 

1 

i!ij 

V 

” h 

3x4, 


2*0 


B, 


B, 


B, 


(18b) 


(18c) 


(ISd) 


8 


(tSe) 



(180 


The above equations may be compared with similar equations in Weinman and Davies (1978) 
and Liou and Ou (1979). The difference lies in the parameters h !tnd h* in the above equations. 
In this form, the solution for horizontally finite clouds tends to the two-stream solution when 
the sides of the cloud are increased to infinity. It should be noted that althou^ the x- and 
y-directions are interchangeable, that is not true of the vertical direction (note the difference 
between h and h'). This is a drawback as will be demonstrated later when numerical examples 
are presented but the errors are small. In our opinion, this is a modest price to pay for avoiding 
the unrealistic results obtained from the Eddington approximation. 
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4. 


MONTE CARLO SIMULATION 


In order to test the accuracy of the two-stream approach, the transfer problem was also 
solved using a direct Monte Carlo simulation. The Monte Carlo model used is very similar to that 
described by Davies (1978) and Weinman and Davies (1978), and only the extensions needed for 
the current application are given here. 

It is helpful to consider tluee separate sources of theimal photons: those emitted uniformly 
from within the doud, those emitted from the ground immediately below the doud base, and 
those emitted from the ground which strike the cloud sides. A summary of the simulation 
characteristics of each source type is, presented in Table 1, which shows the energy per simulated 
photon as well as how the choice of Initial coordinates and direction cosines for each photon is 
made. 

As a result of the subsequent random walk of each photon, simulnting the transfer of ther- 
mal radiation, a number of photons from each source type may ultimately exit the cloud surface 
throudi an area A in a solid angle bin at a mean angle Q to the normal to A. That is 

/*®2 r^i 

cos 8 = I I cos 0 sin dfi d 0 (19) 

where 9i, 62 * ^ limits of the solid angle. If n^ such photons have energy each, 

02 have E 2 , and 03 have E 3 , the total simulated radiance through A in the direction defmed by 
$ is then given by 

I * (Hi El 02 E 2 + n3E3)/Acos0 (20) 

The stochastic nature of the Monte Carlo simulation causes its results to contain an element 
of uncertainty, which may be objectively assessed (Davies. 1978). We also note that the Monte 
Carlo solutions tend to the correct plane parallel solutions as s /2 -► o®. 


As an additional check of the two-stream approximation, the source of photons incident on 
the cloud sides was altered in the so-called “modified Monte Carlo*' version. For this version, 
the energy per unit area incident on the cloud sides was unchanged, but the initial direction 
cosines were chosen to be equally distributed in the upward and downward directions. This 
simulated the two-stream approximation more closely, but it no longer directly simulated the 
actual problem. In the next section, two-stream results will be compared with Monte Carlo 
simulations. 
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5. RESULTS 
a. Fluxes 

The solution of (16) subject to the boundary conditions (18) is obtained using finite Fourier 


transforms and is described in detail in Weinman and Davies (1978). We start by defining 


( 21 ) 


2 2 


where Ci and are eigenvalues to be defmed later. All the quantities of interest may be derived 
from Iq which is given by 


® ° i j IfiS + sin (fiS)! l?js + sin (?jS)l ^ 2 


( 22 ) 


Here, J’j and I’j satisfy 


and 


where 


and 


ri,j * hcot(ri,jS/2) 

^ = C+e^2 + C.e-^z - (x/a2) 
A2 a \2 + f? + f? 


(23) 


(24) 


(25) 


(26) 


The constants C^. and C. are determined from the transformed boundary conditions (1 8a) and 
(18b) which are 


and 


\ h' dz/,.,. rifj V2/ \2/ 

* sm ( — sin — 

V h' 3z/,.o fifj V:/ \2j 


(27a) 


(27b) 
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To obtain the flux exiting each face of the cloud Eqs. (17aH17h) are evaluated at the ap- 
propriate face. The monochromatic emittance, transmittance and reflectance of plane parallel 
clouds have been used to identify their radiative behavior. Analogous deflnitions for horizontally 
finite clouds are not universally accepted and there is no general consensus available on these 
parameters. Here we deflne the emittance in terms of the flux exiting a particular face of the 
model cuboidal cloud. For an isolated, isothermal cloud at temperature Tq, the emittance of the 
top face is the ratio 

s 

and the emittance of the side face at x » ~ is 

2 

and so on for all the faces. Note that there is no ground contribution and F 2 +, j.*. is a function 
of x and y while F^^., x-t* ^ ^ function of y and z, unlike the plane parallel problem. Now if we 
consider the cloud at absolute zero over a black ground at temperature Tp the reflectance of 
the bottom face of the cloud is defined as the ratio 


(Fz-. z-),.o/'8. 

The transmittance of the top face of the cloud is defined as the ratio 

and we have also defined the flux exiting the four faces as transmittance, i.e. 

and so on. *nie advantage of these definitions is that they coincide with the plane parallel 
definitions when the sides of the cloud extend to infinity. Moreover, for a black ground 
the solutions for the cloud at absolute zero and ground at absolute zero can be superim- 
posed to yield the radiative Held. The present definitions are only a function of cloud 
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characteristics, i.e., SS, g, and optical dimensions z* and s. They are thus independent of tempe^ 
ature and they differ from those presented by Liou and Ou (1979). 

The point hemispheric flux exiting a cuboidal cloud of optical dimensions (s, s, z*) (1, 1, 1) 

is shown schematically in Figure 2. Scattering parameters CS • 0.638 and g > 0.86S are represent* 
ative of a C-1 water cloud at lOfim (Deirmendjian, 1969). In the rigures, H and L refer to rela* 
tively high and low flux and the contours indicate the flux pattern. The notable feature here is 
that emission is strongest from the interior regions of the cloud face whereas upward transmission 
is greater around the edges. The source of transmitted radiation is below the cloud so that 
relatively little radiation from below penetrates to the top and there is a gradient in the flux 
exiting the sides. The flux patterns for a superposition of an isothermal cloud on a black 
ground will be the Planck weighted addition of the emission and transmission results. Thus the 
composite case may show edge darkening or brightening of the top face depending on the optical 
dimensions, scattering properties and the temperatures of the cloud and ground. 

The particular case of a cubic cloud of optical dimensions (s, s, z*) (10, 10, 10) at 250K 

above a black ground at 300 K is shown in Figures 3 and 4. Figure 3 gives the upward flux exit* 
ing the top face, (Fj+^ 2 +)^,^*, expressed in equivalent black body temperature for both the iso- 
lated (no ground) and composite configurations. The flux in the central portion of the cloud 
face at optical distances X’* (» 1.43 for this case) from the edge, is equal to the value obtained 
for plane parallel clouds. Figure 4 shows the upward component of the flux exiting each side 
of the cube, “ **ven by (17c). Note that the effect of the underlying 

base can be seen in the composite case at the bottom edge. The figures show quite clearly that 
finite clouds cannot be considered black emitters when the radiative flux is computed. 

A comparison of fluxes obtained with the Monte Carlo method is shown in Figure 5; the 
Monte Carlo results are shown in parenthesis. Values correspond to the average flux in the 
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boxed area expressed in with the lead digit 2 suppressed for darity. The boxes left unfilled may 
be completed using symmetry. The largest differences appear at the edges in the composite 
case (->• lO^K) and may be attributable to the fact that the diffusion approximation breaks 
down. An error is also introduced in the comparison over large bins where the flux varies 
rapidly. The RMS differences between the analytic results and Monte Carlo simulations are 
summarized in Table 2 in which the (10, 10, 2) case has also been included. It may be 
noted that use of the ‘^modified Monte Carlo** reduces the RMS error for side fluxes con- 
siderably. This indicates that there is an inherent defect in the use of flux boundary condi- 
tions in Eqs. (IScHlSO. 


The flux patterns shown averaged over each face yield the mean hemispheric flux exiting 
each face. This quantity may then be used to perform a bulk radiative balance for the cloud 
as a whole. From Eqs. (17aH17h) the average fluxes may be written as 


for the top face. 


2 2 




2 2 


for the bottom face, 


' *- 8 / 2 “ -^0 /i ^ ^ ^’‘/x-s/2 


for one side face, and so on. 
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The results of the integration, using (22) are as follows 


m ^ ^ 16 sin (fj s/«) sin (jj s/2) 

I : 1 j lfi> + »•" (fi»)l ({jS + sin (fjS)l 


( 28 «) 


[‘ 


I 

rsi I 


(28b) 




(28c) 


and so on. 

Thus the mean emittance of the top face of an isothermal, isolated * 0) cuboidal cloud 
is the ratio (F 2 +,z+)j,j*/»Bo. This is plotted in Figure 6 for two finite cloud geometries, 
one cubic and one with horizontal dimensions ten times the vertical; the ordinate shows the 
vertical optical thickness, z*. Also plotted are plane parallel results obtained using the two- 
stream approximation and an accurate Neumann solution using the Henyey-Greenstein phase 
function.^ The mean emittance of the top face of the cuboidal cloud falls short of that of 
plane parallel clouds of the same optical properties and depth. This differs from the results 


S. Olson (porsonai communication) 
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of the hddington approximation shown in Figure 6 of Uou and Ou (1979) in which the (lOr* 
lOr, r) case shows emittance in excess of the plane {^ntlel cloud. Finite 6 of the present paper 
shows that Eddington emittance is consistently larger than the exact results and that it exceeds 
unity when z* ^ 5. Thus the finite Eddington approximation will overestimate the emittanw, 
while the two-stream is seen to be an underestimate, but the plane parallel result are very clme 
to the exact for z* 1. 

Although the top face of the model cloud emits less than the corresponding plwe parallel 
cloud, the total emission from the cloud must include the contribution of the sides of the cloud. 
This is shown in the flgure as the upward component of the side face emittance. For the cubic 
cloud, the contribution of the four side faces exceeds the emission from the top face and the 
total emission from the cloud is actually greater than the corresponding plane parallel cloud. The 
(lOr, lOr, r) case also shows total emittance greater than the plane parallel. If the cloud is not 
isolated, a portion of the side emission will be intercepted by neighboring clouds, thereby de> 
creasing the total emittance. 


b. Cooling Rates 

The net flux divergence within the cloud gives the infrared cooling or heating experienced 
by each element of the cloud. For the three dimensional problem the cooling rate is given by 


3T I 

dt * pCp ^ 


(29) 


where the net flux divergence is obtained from Eq. (II). Here p is the air density assumed to be 
7.2 X I0”*gcm*^ and Cp the specific heat at constant pressure taken to be I J g“‘ ®C'^. For 
illustration, we have computed spectral cooling rates in the 8-13.6 pm region. It was shown by 
Stephens (1978) that the window cooling rate forms the bulk of the total cooling rate for pbne 
parallel clouds, so the results presented here may be considered representative of the total flux 
divergence. The cloud is cubic of dimensions I km on each side at a temperature of -lO^C, the 
ground is at 15*C. It is further assumed that the cloud i? type C-l of droplet density 100 
particles/cm^ having optical properties tabulated in Deirmendjian U969), 
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Figure 7 compares the horizontdiy averaged cooling rate of the cubic cloud with a plane 
parallel cloud of the same height. The additional cooling from the sides of the cloud increases 
the mean cooling at each level in the cloud. The level at which cooling to space is overcome by 
heating from the ground is lowered by 330m. What is even more striking is the distribution of 
net flux divergence at each level. Figure 8 shows the heating and cooling at four levels within 
the cloud marked A, B, C, D on Figure 7. The stippled areas indicate that heating is confined 
to the core at the 400 and 300 m levels, and cooling is around the edges. In regions close to the 
doud base, such as level D, which is 100m above the base, the horizontally aven^d divergence 
indicates heating as is evident from Figure 7 but Flgu**e 8 shows that there is an outer skin that 
exhibits cooling. 

The model illustrated here has a low droplet density. Additional computations have been 
performed for denser clouds and are presented in Table 3 which shows the horizontally averaged 
net flux divergence at selected levels within the cloud. Optically thicker clouds have more ex* 
treme divergence rates at the top and base, consistent with results obtained from the plane 
para" nodel. 

c. Radiances 

Although diffusion approximations furnish respectable results for fluxes,, the radiance field 
derived from such methods is totally inadequate. For example, the two-stream approximation 
assumes a uniform radiance in each hemisphere, however, the two-stream approximation can 
be used to compute the scattering source function. J,, from which the radiance emerging from 
the cloud boundary can be computed by path integration (Weinman and Davies, 1978). 

The radiance exiting a finite cuboidai cloud as seen by a sensor that measures the area 
average radiance at some viewing angle is shown schematically in Figure 9. It is assumed that 
the radiance for azimuthal angle, 0 * 0 . will suffice for comparison with a real cloud of slightly 
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irregular shape. Because it is impossible to represent the detailed features of a real cloud, one 
azimuth angle is as good as any other for the level of sophistication of the model. The radiance 
exiting at a point xo, Xq. z^, on the face of the cloud may be written as 

I(xo.yo.*o*®*0)“Bie'S*(«)+ ((l-a)Bo+J,(xo,yo,2o(x.^).^.OlV"*^®^“ (30) 

/c/Cd) J V 

where J| is the scattering source function given by (3) and the slant optical path lengths S(d) 

and S*(d), the limits of integration x^(d) and x/(d) and the functional dependence z(x, 9) which 

defines the slant path are tabulated in Weinman and Davies (197$) and will not be repeated here. 

Note that the sign convention in that paper is the reverse of that used in this study. Davies 

(1980) has shown that this method of computing exiting radiances is superior to that which a 

first order Eddington expansion provides. 


it remains to define the scattering source function in the two stream approximation. Here 
we may use the complete phase function expansion from Giandrasekhar (i960) instead of (6). 

p(M,d;M'.^')- £ aS/(p/(M)P/(M’) *2 y — ~ Pf”(M) P,"*(m') cos m (d -d')\ (31) 

So V mil m)\ ' ' } 

where 0/ are the cocfncicnis in the phase function expansion and P/(/x) and PJ"(m) are the 

Legendre functions. The scattering source function can then be written as 

X iX '*’X <***'! «*^' (32) 

Now substituting for I* from (8) and p(fi, d; m’. d*) from (31), Eq. (32) can be integrated analyti* 
cally using relations available in Hobson (1955) ano Gradshteyn and Ryzhik (1965). The result is 

^ _ r(-l) (-3) (-5). ..(-:/+ 1)1 'I 

J. . 01, * 01. {— + Z 0:,.. L,,,. 


cSir 

+ — - vdjj cos d + ly sin d) 
8 


(33) 


|o, ♦ J 0;,.. i 
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where 


jFj - I + 


r(2)r(2) 






T(2 + k) r(2 + k) k! 


If we choose the Henyey-Greenstein phase function to represent p, then 23^ ■ 3g and 232/4.1 ■ 
(4/ 3) <i< Hulst, 1980). For zenith angles less than 70* the series can be approxi- 

mated quite well by 


J, • 23 + gl| d + ly sin 


which has been used with a slight modification by Weinman et al. (1981). 


(34) 


After substituting (33) into (30), the radiances obtained can be analytically averaged over 
the projected anas furnishing lengthy expressions for the desired quantity. An illustrative result 
is shown in Figure 10 in which the emergent radiance of an isolated cuboidal cloud is plotted 
against the zenith angle of the viewing sensor. The radiance values are normalized against the 
black body emission of the cloud, Bq. The case z* « S is considered; individual curves are labeled 
with the horizontal optical thickness, s. The anisotropic nature of the emeiient radiation is evi- 
dent. The radiance is sharply reduced for finite clouds and there is a significant angular depend- 
ence. It may be noted that the curve for (5, 5. 5) is nearly symmetrical about 6 * 45®. This 
shows that the slight asymmetry alluded to in Section 3 is of little consequence numerically. 

Also marked on the figure are plane parallel solutions obtained by path integration of tlie Edding- 
ton source function and a 16-stream discrete ordinate solution.^ The two-stream solution is satis- 
factory for all angles while the Eddington approximation yields poorer results e.xcept around 50®. 


■■■■ 

*M.*L. Wu (petsonal commumcition). 
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The above case considered an isolated cloud. The combined radiance profile of a cloud at 
2S0^K overlying a black surface at 300*K is shown in Figure 1 1 for optiol dimensions (10, 10, 
2) and (10, 10, 10) compared with the plane parallel solutions. The rontribution of the side 
faces to the radiance at laiye zenith angles is a prominent feature of the results. The local maKi* 
mum is a function of the geometry of the cloud. It may be noted that even at d ■ 0* the finite 
cubic cloud with z* ■ s ■ 10 appears warmer then the plane parallel model suggests. 

This discrepancy can a^bct the determination of cloud top heights from the observed brightness 
temperature of small cloud elements. The solution here assumes that the ground is at a uniform 
temperature over a large distance compared to the cloud dimensions. It is aho possible to solve 
problems where the ground temperat^xre immediately below the cloud differs from the surround* 
inp, (see Weirunan et al., 1981). 

Another observation from Figure 1 1 is that the two-4tream approximation is good when 
the optical depth is large but that it fails to reproduce the radiance profile with z* * 2 in the 
plane parallel ose. This is not a serious handicap beause clouds in the infrared window usually 
have optical thidcness of -»10 or more. 

The detailed radiance pattern as viewed at 8 ■ 50^, p ■ 0* is shown in Figure 12 in units of 
brightness temperature for a cubic cloud of optical dimensions (10, 10, 10). The pattern is 
shown for both the isolated cloud at 250*K and for the cloud over a black ground at 300” K. 

The sharp gradient at the edges of the cloud are a feature of the results, although the bulk of the 
cloud is seen to exhibit a nearly uniform radiance. Figure 1 3 shows a comparison of the analyti* 
cal results with Monte Carlo simulations (in parenthesis) for the same case. Again, the lead digit 
2 has been suppressed for clarity. Errors are seen to be confmed primarily to the edgwi where 
the Monte Carlo results indicate a sharper gradient than the analytical results. The overall com- 
parison is summarized in Table 2. 


The same model cloud when viewed from directly overhead at 9 ■ 0* shows the radiance 
pattern in Figure 14. It may be noted that the central portion of the cloud is bri^ter in emis- 
sion (isolated) whereas the composite cloud on black ground shows bri^tening about the edges. 
Comparisons with Monte Carlo simulations are presented in Table 2. For the radiance, as for 
flux, the '‘modified Monte Carlo" generally gives closer agreement with the analytical results, 
again highUehtii^ the problem introduced by the side flux boundary conditions of Eqs. (I8c)-(i80> 


6. SUMMARY 


The results presented here illustrate the differences between radiation from finite and plane 
parallel models of clouds. These are 

(1) Finite clouds do not behave as black emitters if their horizontal extent is of the same order 
as their vertical dimensions unless the optical thickness is much greater than the diffusion 
lengrh, 

(2) There is a significant infrared cooling out of the sides of finite clouds. 

(3) The brightness temperature measured from a finite cloud by a low resolution sensor is 
higher than the cloud top temp<*rature. ■ 

(4) The radiance pattern across the face of a fmite cloud shows sharp gradients near the edges. 

The two-stream approximation used in the <tudy was shown to be within a few degrees of 
the Monte Carlo simulations for radiance measurements indicating that this computationally rapid 
method may be advantageously used to model the radiative characteristics of fmite clouds. 
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Table 1 

Monte Carlo Simulation Characteristics 


Energy Pftr Photon 

Starting Position 

Starting Direction Cosines 

(a) Internal Source 

4it(1-S3)BosV/Nio 

uniform within cloud volume 

isotropic 

(b) Source at Cloud Base 


uniform on cloud base 

v sRiW2 


hj * (1 -v2)t/2 cos 2ir Rj 

hj s (1 - sin 2ir Rj 

(c) Source at Cloud Sides 

2ffBi sz’/Njo 

uniform on cloud sides 

hi = Ri^/2 


h 2 » (1 -h?)^^^ cos It R 2 
v » (1 - sin IT R 2 

(d) Modified (c) 

2irBi sz*/N3o 

uniform on cloud sides 

h, « R,''2 


h 2 = (1 cos 2it R 2 

v = (I - h|)‘^2 sin 2ir R 2 


Bo * Planck function at cloud temperature 
Bj *» Planck function at ground temperature 
R 2 are random numbers, uniform on (0, 1) 

V is the vertical direction cosine, positive upwards; h^, hj are the horizontal direction cosines 
Nio, N 20 , N 30 are the total number of simulated photons for each source type 


Table 2 

Comparison Between Two-Stream and Monte Carlo Simulations 
To - 250*K Ti - 300*K 

RMS Difference in 


( 10 . 10 . 10 ) ( 10 . 10 , 2 ) 



Upward 

isolated 

0.7 

0.7 


Hemispheric 

Flux 

composite 

5.8 (MO 

4.8 (MO 


4,9 (MMC) 

3.8 (MMO 


Radiance 

isolated 

2.3 

1.9 

Top 

at * 50* 





n 

composite 

4.5 (MO 
1.7 (MMC) 

2.0 (MO 
2.3 (MMC) 


Radiance 

isolated 

3.6 

2.6 


at a * 0* 

composite 

1.3 (MO 
3.8 (MMC) 

3.4 (MC) 
2.7 (MMC) 


Upward 

Hemispheric 

isolated 

2.0 

7.2 (MO 
3,6 (MMC) 

4.3 

' 5.7 (MC) 
4.7 (MMC) 

Side 

Flux 

composite 



Radiance 
at a « 50* 

isolated 

2.7 

3.3 



H 

composite 

2.1 (MO 

1.2 (MMC) 

2.1 (MC) 
1.8 (MMO 


isolated - cloud alone 

composite - cloud over black ground ‘ 

(MO - correct Monte Carlo version with side flux only from the bottom 

(MMC) - modified Monte Carlo, with same boundary conditions as analytic model, i.e. distributed 
side flux 



Table 3 

Horizontally Averaged Heating (*)') and Cooling (-) 
Rates in ”C/hr in the 8-13.6piin Region 


Height 

From Base (m) 


Droplet Density No./cm^ 


100 

200 

300 

400 

1000 

-5.3 

-10.8 

-16.2 

-21.6 

950 

-3.0 

-3.5 

-3.1 

-2.5 

900 

-1.8 

-1.4 

-0.9 

—0.6 

500 

-0.4 

-0.4 

-0.4 

-0.4 

100 

+0.5 

+0.2 

+0.1 

+0.2 

50 

+1.2 

+1.5 

+1.2 

+0.9 

0 

+2.6 

+5.6 

+8.6 

+11.6 


Goud is type Cl with horizontal dimensions ( I km, 1 km, I km) Tq * -10®C; T j * 1 5®C 






Fig. 1. Schematic ii 
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CUBICCLOUD(10, 10, 10) TOP 
To = 250«K Ti=300®K 


ISOLATED I I COMPOSITE 



Fig. 3. Upward hemispheric flux, expressed as effective black body temper- 
ature, exiting the top face of a cubic cloud of optical dimensions 
(10. 10. 10). Isolated - cloud alone; composite - cloud over black 
ground. 




CUBIC CLOUD (10, 10, 10) SIDE UP 
To = 250»K Ti=300'»K 

ISOLATED ! ! COMPOSITE 

I I 
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Fig. 4 . Same as Figure 3, but upward component of hemispheric Hux exit- 
ing each side face. 
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Comparison between two-stream and Monte Carlo simulations (in 
parenthesis) for cases given in Figures 3 and 4. Flux is expressed 
in effective black body temperature with the l*ad digit Z suppressed 
for clarity. Empty boxes may be filled by symmetry. 
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VERTICAL OPTICAL THICKNESS 

Fig. 6. Mean emittance of cloud faces. Henyey-Greenstein phase function 
has been used for the Neumann solution. 
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COOLING {-) AND HEATING ( + ) RATES "C/hr 


Fig. 7. Horizontally averaged heating (•►) ajid cooling (-) rates in the S 
13.6 Mm region within a 1 km cloud at -10® C overlying a blac 
surface at 15® C. The cloud is type C-1 and has 100 particles, 
cm^. Stippled areas show heating. 



Fig. 8. Horizontal section at 400m (A), 300m (B). 200m (C) and 1 00m 
(D) above the base of the cloud showing the heating (*♦•) and 
cooling (-) rate distribution in C hr at each level. 
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Fig. 9. Viewing geometry- for sensing average radiance from cuboidal clouds. 




NORMALIZED RADIANCE I/Bq 


1.0 



ZENITH ANGLE(DEGREES) 

Fig. 10. .Average normalized radiance as a function of zenith angle for an 
isolated cloud of optical depth, z*=5. Individual curves are 
marked with the horizontal optical thickness, s. .Also shown are 
plane parallel results using the Eddington approximation and the 
sixteen-stream solution for which the Henyey-Greenstein phase 
function was used. 


BRIGHTNESS TEMPERATURE (“ K) 



0 10 20 30 40 50 60 70 80 90 

ZENITH ANGLE (DEGREES) 

Fig. 1 1. Average raJiance expressed as the brightness temperature as a 
function of zenith angle for clouds of optical dimension (10. 

10. 2) and (10. 10. 10) compared with corresponding plane 
parallel clouds. Tq= 250“ K. Tj= 300* K. 
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Fig. 13. Brightness temperature companson between two-stream and 
Monte Carlo simulations (in parenthesis) for case shown in 
Figure 12. Lead digit 2 suppressed for claritv. 














































































CUBIC CLOUD (10, 10, 10) 
BRIGHTNESS TEMPERATURE AT 6 = 0 ° 
To = 250°K 
Ti = 300°K 
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